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The recent detection of a convective core in a main-sequence solar-type star is used here to test 
particular models of dark matter (DM) particles, those with masses and scattering cross sections in 
the range of interest for the DM interpretation of the positive results in several DM direct detection 
experiments. If DM particles do not effectively self-annihilate after accumulating inside low-mass 
stars [e.g. in the asymmetric DM scenario) their conduction provides an efficient mechanism of 
energy transport in the stellar core. For main-sequence stars with masses between 1.1 and 1.3 M©, 
this mechanism may lead to the suppression of the inner convective region expected to be present 
in standard stellar evolution theory. The asteroseismic analysis of the acoustic oscillations of a star 
can prove the presence/absence of such a convective core, as it was demonstrated for the first time 
with the Kepler field main-sequence solar-like pulsator, KIC 2009505. Studying this star we found 
that the asymmetric DM interpretation of the results in the CoGeNT experiment is incompatible 
with the confirmed presence of a small convective core in KIC 2009505. 


I. INTRODUCTION 


The quest to unveil the constituents of the Dark 
Matter (DM) of the Universe is nowadays one of the 
most challenging goals in modern physics [l|. Several 
models beyond the standard model have been put 


forward to explain the properties of DM in terms of new, 
still undetected particles [ 1 , |^. Colliders, ground-based 


telescopes, satellites, and underground experiments are 
trying to detect the hypothetical non-gravitational in¬ 
teractions of these DM particles aiming at the discovery 


of their elusive nature [J-l7| . Interestingly, several direct 
detection experiments have reported signals which could 
have arisen from collisions of low-mass DM particles 
with the nuclei in the detectors (MU. However, this 
interpretation is challenged by the null results in many 
other experiments [im3- Whether these opposed 
results axe fully incompatible or not is still a matter of 
debate 


Here we apply a novel approach to test the existence 
of DM particles with properties as those that would ex¬ 
plain the positive results of direct detection experiments 
by looking for the particular effects that they would 
produce in stars slightly more massive than the Sun. 
As it will be shown, this method offers an inexpensive, 
complementary probe on the plausibility of some DM 
candidates. The models of DM that can be tested with 
our approach are non-(or feebly) annihilating particles, 
with the low masses and relatively large scattering 
cross-sections with baryons required to explain the 
results of the DAMA Q, CoGeNT @ and CDMS-Si [ll| 
experiments. Theories that predict these low-mass 
Weakly Interacting Massive Particles (WIMPs) include 
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Asymmetric DM (ADM, [2ll - [^ ). although our analysis 
applies to any model of DM for which the number of 
DM particles accumulated in the core of a star is not 
significantly depleted due to annihilations or decays, 
independently of the mechanism behind this property. 
In the context of ADM, although the self-annihilation 
cross-section of the DM particles may be large, their 
annihilation inside stars is prevented by the lack of 
antiparticles, originated from an asymmetry in the par¬ 
ticle/antiparticle abundances related to the baryonic one. 

The properties of ADM have already been constrained 
using observations of other type of stars. Neutron stars, 
being very dense, could capture and accumulate huge 
quantities of DM particles ^ , leading to changes in 
the properties of the star [ 2 g-[^ and the hypothetical 
creation of a black hole that could destroy it [29l - [3^ . 
Although the constraints on ADM from this approach 
are very stringent, they have the uncertainty associated 
with the complexity of the high-energy physics involved 
in the self-collapse of DM in the core of a neutron 
star [^. 1 ^. 


On the other hand, most of the processes governing 
the interior of main-sequence stars are well understood 
and tested by precise observations, including the solar 
neutrinos and helioseismology [4l|. In the Sun, the 
accumulated ADM particles would quickly thermalize 
in a small region of the solar core and their conduction 
would modify the solar central temperature and density, 
leaving potentially stro ng i mp rints in the solar neutrino 
fluxes and oscillations |4 ^h 4^ . Remarkably, stars with 
masses slightly greater than the Sun can have stronger 
structural changes in their interior due to ADM. It has 
been shown that if DM does not self-annihilate after 
accumulating inside stars with masses between 1.1 and 
1.3 M 0 , their conduction may remove energy from the 
stellar nucleus efficiently enough to suppress the small 
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TeS log(9') [^/^Is -Rcc (Ai/)oi 2 |S'{At/o yoio}! 

LiU( R* ) ( mHz ) 

- Observations: 

Dushera 6200 ± 200 4.30 ± 0.2 0.019 ± 0.006"- 0.061-0.071*’ 88 ± 0.6 0.0032 ± 0.0006 


- Models: 

Standard (no DM) 6222 ± 107 4.21 ± 0.01 0.019 ± 0.004 0.069 ± 0.006 87.3 ± 1.1 0.0028 ± 0.0011 
DM model 1 6208 ± 110 4.21 ± 0.01 0.019 ± 0.004 0.069 ± 0.006 87.3 ± 0.9 0.0028 ± 0.0012 

DM model 2 6208 ± 111 4.21 ± 0.01 0.019 ± 0.004 0.069 ± 0.006 87.2 ± 1.1 0.0028 ± 0.0012 

DM model 3 6255 ± 103 4.21 ± 0.01 0.020 ± 0.004 0.001 ± 0.008 87.1 ± 1.0 0.0006 ± 0.0004 

DM model 4 6219 ± 107 4.21 ± 0.01 0.020 ± 0.004 0.028 ± 0.016 87.3 ± 0.9 0.0028 ± 0.0015 


^ the observed [Fe/H]= log(Z/X) — log(Z/X)0 was converted to the surface \Z/X] assuming equal iron and metal abundances, 
the range in the convective core radius is that of the best fit models in Ref. 

TABLE I. Parameters determined from ground-based high-resolution spectroscopy and Kepler high-precision photome¬ 
try for Dushera, and distribution of the same parameters in the models obtained from our grids. 


convective core that these stars are expected to develop 
in the main sequence (4^ . 

In the last years, the space-based asteroseismic 
missions CoRoT and Kepler have detected solar-like 
oscillations in thousands of stars [l^ . The high precision 
in the determination of the oscillation frequencies of 
these stars has allowed not only measurements of their 
masses and radius with unprecedented accuracy, but also 
opened a window into the stellar interiors, including the 
discrimination between He-burning and inert cores in red 
giants [i^. In this work, we analyse KIC 12009504, the 
first Kepler main-sequence solar-like pulsator to have its 
convective core detected thanks to asteroseismic studies 
performed on this star [s^. We use this detection to 
rule out DM models that would prevent the formation 
of such a convective region. 

This article is organized as follows. The modeling of 
the star and the DM particles are described in Section HIl 
Section nil reviews the impact of DM conduction in 
stars. In Section |TV] we present the diagnostic tools used 
to infer the presence/absence of a convective core from 
the stellar acoustic oscillations. In Section El we show 
the results of our simulations and how they compare 
with the observations of KIC 12009504, and finally in 
Section ED we summarize and discuss the conclusions of 
this work. 


II. MODELING 
A. Stellar modeling 

KIC 12009504, also known as Dushera, is a main- 
sequence solar-like pulsator which has been observed 
by ground-based high-resolution spectroscopy and 
by the high-precision photometry of the Kepler mission 


[40| . These observations led to measurements of its 
effective temperature Teff, surface gravity log(( 7 ) and 
surface metallicity [Z/Ajs which are shown in Table |D 
For the modeling of Dushera we assumed the large 
uncertainties in Teff quoted in Table [D following the 
approach in a previous thorough analysis of this star [s^ , 
due to an inconsistency between different spectroscopic 
analysis (d^ . 

Solar-like oscillations have been observed in Dushera, 
with the individual frequencies of 34 modes precisely 
determined in Ref. (50j |. allowing an accurate estimation 
of its properties. In particular, its mass has been 
estimated with a very small uncertainty (±0.03 Mq) 
through the htting of individual frequencies on the 
basis of stellar models [5l|. However, for our model 
fitting we adopted a conservative approach and con¬ 
sidered the broader mass range obtained by Ref. [3^ : 
1.23 ± O.12M0. The latter uncertainty comes from a 
combined analysis of the modeling of Dushera by several 
teams, using different stellar evolution codes, pulsation 
codes, and htting techniques, including the htting to 
frequency ratios which are sensitive to the stellar interior. 

The modeling of Dushera was performed using a 
modihed version of CESAM stellar evolution code 
including microscopic diffusion [5^, convection through 
the mixing-length theory [53| with overshooting of the 
convective core implemented as dov = ctoi; niin(i?cc, 
where dov is the overshooting distance, i?cc is the radius 
of the Schwarzschild convective core, H„ the pressure 
scale height and aov a free parameter Im, and stellar 
heavy element mixture as the solar one [56j |. 

An extensive grid of stellar models was computed, 
with stellar masses in the range M* = 1.1 — 1.35 Mq, 
initial metallicities between Z = 0.012 — 0.024, initial 
helium abundances between Yq = 0.26 — 0.30, over¬ 
shooting with efficiency parameter aov = 0.05,0.15, and 







3 


convection with aMLT =1-6, 1.8 and 2.0. The same 
grid was computed for the standard stellar evolution 
scenario, without DM, and also considering the impact 
of the different models of DM particles described in 
Table in and in Section Pi B1 

From all the stellar models computed in our grids, we 
selected those that reproduce the observed properties 
of Dushera (Teff, log((?), \ZlX\s) within l-cr and the 
large frequency separation, Az/, within 2-cr. We allowed 
for a larger error in ISv because, at this stage of the 
modeling pipeline, this parameter was estimated using 
scaling relations based on an asymptotic approximation: 
^ (see 

[ 5 ^). Our grid-based selection procedure is the 
standard to determine stellar properties from their 
oscillations (see e.g. Ref. [s^) and is comparable to 7 
of the 8 methods compiled in Ref. [s^ to simulate the 
same star. An alternative would be to use Markov Chain 
Monte Carlo or genetic algorithms, which can be more 
precise in providing the best-fit model due to the finite 
resolution of grid-based methods [HI, 11^. 

In addition, we estimate that the fraction of valid 
models of Dushera that were discarded due to the 
limits imposed in our grids is below 5% (the mass and 
metallicities of the selected models have mean and 
standard deviations of l.I9±0.05 Mq and 0.0I8±0.003, 
respectively). 

This procedure resulted in more than 30,000 valid 
models of Dushera. For all these models, the acoustic 
oscillation frequencies were computed using the ADIPLS 
code [Hi- With these frequencies we recalculated the 
large separation Av and computed the seismic parame¬ 
ters described in Section HVl 


B. Dark Matter modeling 

The present status of direct dark matter detection ex¬ 
periments is characterized by the tension between the 
promising positive results in some experiments id and 
the robust incompatible limits set by the null results in 
others [iM3- While some models try to simultane¬ 
ously explain all the results, other analysis suggest that 
backgrounds or biased analysis may explain the contro¬ 
versy [ild- In this context, the method we propose 
here provides an alternative test of the existence of the 
low-mass WIMPs with “large” scattering cross-section on 
nucleons that best fit the signals in some of these exper¬ 
iments. 

The properties of the models of DM tested here, 
shown in Table im are those that can explain the 
signals in DAMA, CDMS-Si and CoGeNT in terms of 
spin-dependent (SD) WIMP interactions on protons [l^ 



( GeV ) ( 

;cmM 

Experiment 

DM model 1 

7 

10-40 

DAMA+CDMS+CoGeNTH 

DM model 2 

13 

10-36 

dam43 

DM model 3 

8 

10-33 

CoGeNT^’ 

DM model 4 

10 

10-32 

CDMS-Si'^ 


^ interpreted in terms of SI WIMP-nucleon interactions 
^ interpreted in terms of SD WIMP-proton interactions 11^ 


TABLE 11. Characteristics of the models of dark matter par¬ 
ticles tested in this work. 

(labeled “DM model 2, 3 and 4”, respectively), and 
the model that best fits all these experiments in terms 
of spin-independent (SI) interactions on nucleons [Hl | 
(labeled “DM model I”). The DM scattering off nuclei 
through SD (axial-vector) and SI (scalar) interactions 
is implemented in our models separately, and we tested 
models with either pure SD interactions (DM models 2, 
3 and 4) or pure SI interactions (DM model I). How¬ 
ever, more realistic models may interact simultaneously 
through the two classes of couplings: SD interactions 
with hydrogen nuclei and SI interactions with all the 
elements of the stellar plasma. As stated before, we 
stress that our results apply in general to any model that 
predicts DM particles with masses and scattering cross 
sections similar to those listed in Table HIl as long as the 
number of DM particles concentrated inside stars is not 
notably reduced by mechanisms such as annihilations or 
decays. 

We have analyzed the standard scenario of elastic 
scattering between DM particles and nucleons. For 
the sake of simplicity, here we have not considered 
other scenarios with more complex interactions, such 
as long-range interactions, velocity or momentum- 
dependencies, isospin violation, exothermic or inelastic 

scatterings, etc.j _ although they may boost the DM 

impact on stars [6a-l66j|. In principle, the method of DM 
search proposed here could be applied to constrain the 
parameters describing the WIMP-nucleus interactions 
in more general frameworks such as DM effective field 



The inclusion of the impact of DM on the stellar 
properties was computed following the prescriptions 
of Gould 0 for the capture rate (using modified 
subroutines of the DarkSUSY code [Z^), and Gould and 
Raffelt [U for the energy transport by DM conduction. 
The reader is addressed to the references above for a 
thorough description of the formalisms used to model 
the impact of DM on stars. Nonetheless, it is worth 
highlighting some simple relations that illustrate the 
main dependencies of the processes involved. 

The efficiency of the capture of DM particles by 
a star is proportional to the DM-nucleon scattering 
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cross-section and the density of DM around the star, 
oc cFxPx^ depends on the velocity distribution 

of the DM particles. We assumed a canonical DM halo 
model with a Maxwell-Boltzmann velocity distribution, 
a dispersion = 270 km s“^ and a local DM density 
around the star = 0.4 GeV cm“^, as estimated in the 
solar neighborhood [^. Departures from the standard 
Maxwellian velocity distribution are well motivated by 
simulations and observations and have been 

extensively explored to assess the event rate in DM 
direct detection experiments [t^, [ill- We address the 
reader to Ref. [t^, where the impact in the capture rate 
Cx of the uncertainties in Px and was studied. 

The captured DM particles quickly thermalize with 
the stellar plasma and concentrate in the core of the star 
with a distribution characterized by a radius that shrinks 
with the DM mass: Tx oc Consequently, DM 

particles with low masses can influence a broad region 
transporting energy by conduction (r^ Ri 5% R,^ for 
rux = 10 GeV), whereas DM particles with high masses 
distribute too concentratedly to decisively impact the 
star by this mechanism. 


WIMPs remove energy from the stellar core and bring 
it to a distance characterized by their mean free path in¬ 
side the star lx{f)- This parameter defines, together with 
Tx, the Knudsen number of the system K = lx(0)/rx, 
which shows whether WIMPs transport energy locally or 
not. The luminosity transported by DM is: 


Tx.trans(c) = 47rr^ n^.LTEl^’) lx 


X 


ksTir) 

mx 


dT 

kB^f{K)h{r), ( 1 ) 


where K(r) is the thermal conductivity for a gas of 
WIMPs and the elements of the stellar plasma, found 
by in terp olating in the table calculated numerically by 
Ref. [8^. h{r) and f{K) are suppression factors, the 
latter introduced to extend the expression for the energy 
transport in local thermal equilibrium (LTE) to regimes 
with a larger WIMP mean free path within the stellar 
plasma. The real WIMP distribution inside the star, 
n^(r), is found by interpolating between the LTE and 
the isothermal distributions (see Ref. [Slj for the details). 


III. DARK MATTER IMPACT ON STELLAR 
PROPERTIES 

The additional mechanism of stellar energy transport 
by DM conduction described above was implemented in 
the stellar evolution code CESAM. Our modified version 
of the code has been tested in the case of the Sun and 
other stars [13, Ull producing results in agreement 
with those obtained independently by other groups using 
different codes [Slj, . 


The main impact of DM cooling in the properties of 
main-sequence stars occurs in their core, and consists 
in a reduction in the central temperature of the star 
and an increase in their central density (see Ref. 0). 
Strong deviations in the global properties of the star, 
such as its luminosity and its effective temperature, are 
only expected if the star is located in environments with 
DM densities orders of magnitude higher than those 
considered here [s^. 

In this work we are especially interested in the changes 
on the temperature profile in the core of the star. For a 
star like Dushera, the DM conduction will cool the cen¬ 
tral ^5% of the star, leading to an isothermal core which 
is only present when the impact of DM is included (see 
Figure Ilia)). As it was shown in Ref. this cool¬ 
ing may lead to a strong structural change in the core 
of the star, suppressing convection in the central region. 
Convective instabilities arise when the temperature gra¬ 
dient is so steep that rising bubbles of stellar plasma 
do not cool rapidly enough through adiabatic expansion 
and they continue to rise. Due to the large energy flux 
resulting from the CNO hydrogen burning cycle, main- 
sequence stars with masses slightly greater than the Sun 
are known to develop a small convective core in the stan¬ 
dard scenario of stellar evolution. In contrast, this con¬ 
vective core would be suppressed if DM particles as those 
studied here contribute to cool the stellar core. 

Within the context of the ADM models considered 
here, the DM particles accumulate inside the star 
during its evolution without annihilating. This is a 
crucial factor in order to have enough DM particles to 
noticeably cool the stellar core due to DM conduction. 
Depending on the unknown DM properties, mainly nix 
and (Jx-, the reduction of the temperature gradient in the 
nucleus can be efficient enough to prevent the convective 
instability in the core of Dushera. 


IV. ASTEROSEISMIC SIGNATURE OF 
CONVECTIVE CORES 

The presence of a convective core in the center of a 
star has a strong impact on the propagation of acoustic 
waves. The mixing in the stellar plasma induced by 
convection homogenizes the chemical composition in the 
convective regions, thus creating a discontinuity in the 
border with the radiative regions. This is depicted in 
Figure[T]b), where the radial profiles of the H abundance, 
A(r), of two models of Dushera are plotted. An abrupt 
change in the H abundance can be seen in the DM-free 
model, with a convective core at Rcc — 0.07 R*, whereas 
the changes in X(r) for the DM-influenced model, with 
a radiative core, are much smoother. 


Therefore, a convective core leaves a signature in the 
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0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14 0.16 
r (R* ) 

FIG. 1. a) Temperature, b) hydrogen abundance in mass frac¬ 
tion, and c) adiabatic sound-speed profiles for representative 
models of Dushera with a convective core at i?cc = 0.07 R* 
(black continuous line, corresponding to a model not influ¬ 
enced by DM) and without it (blue dashed line, corresponding 
to a model strongly impacted by DM conduction). 


able, specially from space-based data, other diagnostic 
tools are preferable, such as rgio [^ : 


roio = {roi{n),rio{n),roi{n -b l),rio(n -b 1),...} , ( 2 ) 


where: 


roi(n) = 


rio{n) = 


(r'n-1,0 - 4:iyn-lA + + J^n+l,o) 

8 — t'n-l.l) 

(3) 

-{Vn-1,1 - + GVn.l - dt^n-H.O + ^^n-|-l,l) 

8 {l^n+1,0 — l^n,o) 


(4) 

The frequency derivative of roio may be used to diagnose 
the presence or the absence of a convective core in the 
star and, when present, it may in principle provide infor¬ 
mation about convective core’s properties | 88 l| . 


In this work we considered the slope of the seismic 
parameter (Az/)roio, where Avnj = Vn+i,! — Vn,i is the 
large frequency separation defined as the difference in 
frequency between modes of the same degree Z, and con¬ 
secutive radial order n. It was shown that this frequency 
derivative is sensitive to the frequency perturbation 
induced by the discontinuity in sound speed at the edge 
of a small convective core and can be used to infer the 
amplitude of the relative sound-speed variation at that 
region [13, . 


The absolute value of the slope of the above- 
mentioned diagnostic tool, |iS'{A^roio}|, was computed 
for all models in our grids in the range of the observed 
frequencies [13 using the definitions from Ref. [13 . The 
results of our computations confirm that this slope is a 
good tracer of the presence/absence of a convective core 
in the stellar models. 


radial profile of the mean molecular weight and in the 
pressure modes oscillation frequencies of a pulsating star 
(see [lEHSl)- This signature builds as evolution proceeds 
on the main sequence and is caused by the discontinuity 
in the composition and hence in the sound speed at the 
edge of the growing convective core. This effect is clearly 
illustrated in Figure [TJc), where the sound-speed profiles 
in the inner regions of two representative models with 
and without a convective core are shown. The model 
with a convective core presents a distinct sound-speed 
discontinuity at the border between the convective 
region and the radiative one. 

Particular combination of oscillation frequencies, the 
so-called diagnostic tools, may be able to isolate the sig¬ 
nature left on the frequencies by the presence of a small 
convective core. One such diagnostic tool is dro 2 i 3 , which 
combines frequencies with modes of degree Z = 0 to 3 
[H, 113 • Since Z = 3 modes may not be always observ- 


V. ANALYSIS AND RESULTS 

Here we present the analysis of the properties of all 
the stellar models, computed and selected as described 
in Section m their oscillations, and their comparison 
with the observations of Dushera. The reader can have 
a quick picture of our results by reading Table HI where 
the average and standard deviations of the classical and 
seismic properties of the models are shown. 

In addition, the distribution of the radius of the con¬ 
vective core i?cc and the seismic parameter |S'{Az/roio}| 
obtained in all the models (DM-free and DM-influenced) 
can be easily compared with the i?cc derived in Ref. |38l | 
and the observed |S'{Az/roio}| in Figures [5] and O 
respectively. These figures depict violin plots, a con¬ 
venient method of visually showing the distribution 
of data. The width of the violins is proportional to 
the number of models obtained with a given i?cc or 
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FIG. 2. Distribution of the radii of the convective core, Rcc, in all the computed models that reproduce the observations of the 
star Dushera, in the standard scenario without DM (grey) and including the impact of different DM models (colors, see text). 
The shapes of the distributions (violin plots) show the probability density, and the internal discontinuous lines correspond 
to the first quartile, the median and the third quartile. The yellow area shows the range of the Rcc inferred for Dushera in 
Ref. M. 




xfosCfl) xfz/X], xfsfA^r-ninll 

Xt 

P 

Standard 

0.00 

0.17 

0.00 

0.51 

0.68 

0.88 

DM model 1 

0.00 

0.17 

0.00 

0.53 

0.70 

0.87 

DM model 2 

0.00 

0.17 

0.00 

0.52 

0.69 

0.88 

DM model 3 

0.07 

0.17 

0.03 

18.35 

18.62 

0.0003 

DM model 4 

0.01 

0.17 

0.01 

0.47 

0.66 

0.88 


TABLE III. x^’s and p values of the models of Dushera. 


I^IAr/roioll, whereas the internal discontinuous lines 
show the first quartile, the median and the third quartile. 


For each class of DM model, the x^’s of the properties 
of the stellar models were calculated: 


2 (Qobs Qmod 

Xq = - 

\ ^Qobs 


(5) 


for Q = {reff,log(g), |S'{Az/roio}|}. The in¬ 

dividual Xqj the combined Xt — SqXq: ^-^d the p 
values computed assuming Gaussian errors, are shown 
in Table Hill 


First, we found that the modeling of Dushera without 
DM reproduced with an excellent agreement the results 
of previous analysis on the same star [s^, confirming 
the robustness of our grid and selection procedures. All 
the DM-free models of Dushera were found to have a 
small convective core, with a mean radius of 0.069 R* 
and a standard deviation of 0.006 R*. The seismic 
parameter l^jAi/roio}!, used here as a diagnostic tool 
of the stellar interior, was also found to be in agreement 


with the value calculated from the observed oscillation 
frequencies. 

When the influence of DM particles with the proper¬ 
ties of DM models 1 and 2 were considered, we obtained 
essentially the same results as in the standard scenario, 
i.e the model without DM (see Tables HI and IIIII) . This 
result implies that these models of DM particles do not 
significantly impact the internal structure of the star. 
Accordingly, the distributions of the obtained radii of 
the convective cores found in all the models, as well 
as their seismic properties, are very similar to those 
obtained in the DM-free scenario, as shown in the violin 
plots of Figures [5] and El 

On the other hand, when the influence of the existence 
of DM particles with the properties of model 3 was taken 
into account, we found that only 2% of the models of 
Dushera had a convective core. This is a consequence of 
the efficient energy transport by DM conduction, which 
reduces the temperature gradient in the core of the star 
and prevents convection. 

As expected, this significant modification in the 
structure of the stellar core had a clear impact in the 
stellar oscillations. The diagnostic tool l^jAi/roio}! was 
found to be very sensitive to the absence of a convective 
core in the stellar models. As shown in Tables ID m 
and in Figure^ the parameter |S'{Ai/ roio}| for the DM 
model 3 is clearly in disagreement with the observations. 
Whereas the DM-free and DM models 1 and 2 reproduce 
|S'{Ai/ roioll within one standard deviation, the mean 
|S'{Az/roio}| for DM model 3 shows a discrepancy higher 
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FIG. 3. Distribution of the asteroseismic parameter |S{Ai/roio}| in all the computed models that reproduce the observations 
of the star Dushera, in the standard scenario without DM (grey) and including the impact of different DM models (colors, 
see text). The shapes of the distributions (violin plots) show the probability density, and the internal discontinuous lines 
correspond to the first quartile, the median and the third quartile. The yellow area shows the value of |5'{Ai/roio}| calculated 
from the observed frequencies with its uncertainty. 


than 4-cr raw}\ ^ Assuming Gaussian errors, 

the p values from Xt 8 -'^® P — 0-0003 for this model of 
DM and p ~ 0.88 for the other models, indicating a good 
fit to the data for all the models except the DM model 
3, which can be excluded with a 99% confidence level. 

In the case of DM model 4, the WIMP-proton scat¬ 
tering cross-section is so large that the DM mean free 
path inside the star becomes very small and DM is in 
local thermal equilibrium with the stellar plasma. In 
this scenario the impact of DM conduction is drastically 
reduced, being noticed only by a small contraction on 
average of the convective core size when compared with 
the standard scenario. This impact is not sufficient 
to leave a clear signature in the diagnostic parameter 
|5'{Ai^roio}| and does not lead to a goodness of fit 
significantly different from the DM-free scenario. 


VI. DISCUSSION AND CONCLUSIONS 

Our results show that the existence of asymmetric DM 
particles with the properties that explain the positive 
results in the CoGeNT experiment would lead to the 
suppression of the convective core recently detected in 
the main-sequence star Dushera. We have demonstrated 
that the stellar oscillations, in particular the diagnostic 
tool are sensitive to the stellar core and 

can be used to exclude with a 99% confidence level 
the existence of this particular model of DM particles. 
Whereas the observed non-seismic properties of Dushera 
(Teff, \og{g), [Z/X\s) can be accurately reproduced 


by the DM-impacted stellar models, the parameter 
|5'{Ai/roio}| shows a mean discrepancy with the mea¬ 
sured stellar oscillations above 4 -ct. 

The other models of ADM particles studied here, 
namely those that would explain the signals in DAMA 
and CDMS-Si experiments, were found not to lead to 
significant modihcations in the stellar structure, and 
consequently their existence cannot be disproved using 
our approach. Even though we used the positive signals 
in direct detection experiments to set our benchmark 
models, we note that our approach could be applied to 
any model of DM particles with similar properties. 

In our analysis we included the systematics that arise 
due to the variations of the major standard mechanisms 
that are known to impact the presence or absence of a 
convective core in a star, namely its mass, age, metallic- 
ity, initial helium abundance, and the parameters that 
control the efficiency of overshooting and convective 
mixing. In addition, to grasp the impact of changes 
in the stellar composition, we re-computed our grid in 
the DM-free scenario u sing the old estimation of the 
solar chemical mixture [9l|, resulting in a variation of 
the mean R^c of only 10“^ R*. These results confirm 
the robustness of our approach. Nonetheless, we note 
that unexplored input physics in the modelling may 
introduce further systematic errors (see e.g. Ref [^). 
Regarding the modelling of the DM impact, the main 
additional uncertainty comes from the DM density 
around the star, which we assumed equal to the local 
solar value. The uncertainty in this quantity (see e.g. 
Refs. [H, m) directly influences the DM-proton scat¬ 
tering cross-sections assumed for the studied DM models. 
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